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ABSTRACT 

We investigate the secular evolution of the orbital semi-major axis and eccentricity due to mass 
transfer in eccentric binaries, assuming conservation of total system mass and orbital angular momen- 
tum. Assuming a delta function mass transfer rate centered at periastron, we find rates of secular 
change of the orbital semi-major axis and eccentricity which are linearly proportional to the magni- 
tude of the mass transfer rate at periastron. The rates can be positive as well as negative, so that the 
semi-major axis and eccentricity can increase as well as decrease in time. Adopting a delta-function 
mass-transfer rate of 10~^MQyr~^ at periastron yields orbital evolution timescales ranging from a 
few Myr to a Hubble time or more, depending on the binary mass ratio and orbital eccentricity. Com- 
parison with orbital evolution timescales due to dissipative tides furthermore shows that tides cannot, 
in all cases, circularize the orbit rapidly enough to justify the often adopted assumption of instanta- 
neous circularization at the onset of mass transfer. The formalism presented can be incorporated in 
binary evolution and population synthesis codes to create a self-consistent treatment of mass transfer 
in eccentric binaries. 

Subject headings: Celestial mechanics. Stars: Binaries: Close, Stars: Mass Loss 



1. INTRODUCTION 

Mass transfer between components of close binaries is 
a common evolutionary phase for many astrophysically 
interesting binary systems. Indeed, mass ejection and/or 
accretion is responsible for many of the most recogniz- 
able phenomena associated with close binaries, such as 
persistent or transient X-ray emission, neutron star spin- 
up, and orbital contraction or expansion. Theoretical 
considerations of these and other associated phenomena 
in the literature probe these systems quite effectively, 
yet they often do not consider the effects of any eccen- 
tricity associated with the binary orbit. This can be of 
particular importance for binaries containing a neutron 
star or a black hole, where mass loss and natal kicks 
occurring during compact object formation may induce 
a significant eccentricity to t he bin ary (e.g. Hills 1983; 
iBrandt fc Podsiadlowskill995HKalogerall9 96). After the 
formation of the compact object, tides tend to circularize 
the orbit on a timescale which strongly depends on the 
ratio of the radius of the compact object's companion 
to the orbital semi-major axis. Because of this, orbits 
are usually assumed to circularize instantaneously when 
a binary approaches or begins a mass transfer phase. 

Despite our generally well developed understanding of 
tidal interactions in close binaries, quantitative uncer- 
tainties in tidal dissipation mechanisms propagate into 
the de termination of circul arization timescales. For ex- 
ample, iMeibom fc Mathieul ([2005.) have shown that cur- 
rent theories of tidal circularization cannot explain ob- 
served degrees of circularization of solar-type binaries in 
open clusters. Circularization of high-mass binaries, on 
the other hand, is currently thought to be driven predom- 
inantly by resonances between dynamic tides and free 
oscillation modes, but initial conditions play an impor- 
tant role and an extensive computational survey of rel- 
evant parts of the initial parameter space has yet to be 
undertaken (jWitte fc Savoniidll999ll2001l:IWillems et all 



[200l . 

Furthermore, assumptions of instantaneous circular- 
ization immediately before or at the onset of mass trans- 
fer are in clear contrast with observations of eccentric 
mass transferring systems. In the most recent catalog of 
eccentri c binaries with known a psidal-motion rates com- 
piled bv lPetrova fc Orlovl ()1999r ). 26 out of the 128 Hsted 
systems are semi-detached or contact binaries. Among 
these mass-transferring systems, 9 have measured eccen- 
tricities greater than 0.1. In addition, many high- mass X- 
ray bina ries are known to have con siderable orbital eccen- 
tricities fRagu zova fc Popovll2005 l. While mass transfer 
in these systems is generally thought to be driven by 
the stellar wind of a massive O- or B-star, it has been 
suggested that some of them may also be subjected to 
atmospheric Roche-lobe overflo w at each perias tron pas- 
s age of th e mass iv e donor (e.g . lPettersonlll978|) 



iHuanei (1 19561 ). iKruszewskil ()1964D . and iPiotrowskl 
(pel " were the first to study the effects of mass transfer 
on the orbital elements of eccentric binaries. However, 
their treatment was restricted to perturbations of the or- 
bital motion caused by the varia ble component masses. 
Matese fc Whitmircl (l983l . I1984D extended these early 
pioneering studies to include the effects of linear momen- 
tum transport from one star to the other, as well as any 
other possible perturbations caused by the mass trans- 
fer stream in the system. However, these authors derived 
the equations governing the motion of the binary compo- 
nents with respect to a reference frame with origin at the 
mass center of the binary, which is not an inertial frame. 
Their equations therefore do not account for the accel- 
erations of the binary mass center caused by the mass 
transfer (see i; 13. 3p . 

More recent work on mass transfer in eccentric binaries 
has mainly focused on smoothed particle hydrodynamics 
calculations of the mass transfer stream over the course of 
a few orbits, without any consideration of the long-term 
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evolu tion of the binary l|Lavton et al.lll998t iRegos et al.1 
[20n5h . 

Hence, there is ample observational and theoreti- 
cal motivation to revisit the study of eccentric mass- 
transferring binaries. In this paper, our aim is to derive 
the equations governing the evolution of the orbital semi- 
major axis and eccentricity in eccentric mass-transferring 
binaries, assuming conservation of total system mass and 
orbital angular momentum. In a subsequent paper, we 
will incorporate the effects of mass and orbital angular 
momentum losses from the system. 

Our analysis is ba sed on the seminal work of 
iHadiidemetrioul (|1969b[ ) who was the first to derive the 
equations of motion of the components of eccentric mass- 
transferring binaries while properly accounting for the 
effects of the variable component masses on the stars' 
mutual gravitational attraction, the transport of linear 
momentum from one star to the other, the accelerations 
of the binary mass center due to the redistribution of 
mass in the system, and the perturbations of the orbital 
motion caused by the mass-tra nsfer stream. While the 
equations of motion derived bv IHadiidemetrioul (|1969b[ ) 
are valid for orbits of arbitrary eccentricity, the author 
restricted the derivation of the equations governing the 
evolution of the semi-major axis and eccentricity to or- 
bits with small initial eccentricities. 

The paper is organized as follows. In § 2 and § 3 we 
present the basic assumptions relevant to the investiga- 
tion and derive the equations governing the motion of 
the components of an eccentric mass-transferring binary 
under the assumption of conservative mass transfer. The 
associated equations governing the rates of change of the 
semi-major axis and the orbital eccentricity are derived 
in § 4, while numerical results for the timescales of or- 
bital evolution due to mass transfer as a function of the 
initial binary mass ratio and orbital eccentricity are pre- 
sented in § 5. For comparison, timescales of orbital evo- 
lution due to dissipative tidal interactions between the 
binary components are presented in § 6. § 7 is devoted to 
a summary of our main results and a discussion of future 
work. In the appendices, lastly, we derive an equation 
for the position of the inner Lagrangian point in eccen- 
tric binaries with non-synchronously rotating component 
stars (Appendix and present an alternative deriva- 
tion for the equations governing the secular evolution 
of the orbital semi-major axis and eccentricity assuming 
instantaneous mass transfer between two point masses 
(Appendix |B|. 

2. BASIC ASSUMPTIONS 

We consider a binary system consisting of two stars in 
an eccentric orbit with period Porb, semi- major axis a, 
and eccentricity e. We let the component stars rotate 
with angular velocities fli and ^2 parallel to the orbital 
angular velocity rioib, and assume the rotation rates to 
be uniform throughout the stars. We also note that the 
magnitude of fiorb varies periodically in time for eccentric 
binaries, but its direction remains fixed in space. Because 
of this, the stars cannot be synchronized with the orbital 
motion at all times. 

At some time t, one of the stars is assumed to fill its 
Roche lobe and begins transferring mass to its compan- 
ion through the inner Lagrangian point Li. We assume 



this point to lie on the line connecting the mass centers 
of the stars, even though non-synchronous rotation may 
cause it to oscillate in the direction perpendicular to the 
orbital plane with an amplitude proportional to the de- 
gree of asynchronism (|Matese fc Whitmird 1198 31. Since 
the donor's rotation axis is assumed to be parallel to the 
orbital angular velocity, we can safely assume that the 
transferred mass remains confined to the orbital plane. 

We furthermore assume that all mass lost from the 
donor is accreted by its companion, and that any orbital 
angular momentum transported by the transferred mass 
is immediately returned to the orbit. The mass transfer 
thus conserves both the total system mass and the orbital 
angular momentum. 

We also neglect any perturbations to the orbital motion 
other than those due to mass transfer. At the lowest 
order of approximation, these additional perturbations 
(e.g., due to tides, magnetic breaking, or gravitational 
radiation) are decoupled from those due to mass transfer, 
and can thus simply be added to obtain the total rates 
of secular change of the orbital elements. 

3. EQUATIONS OF MOTION 
3.1. Absolute Motion of the Binary Components 

Following iHadi idemetrioul ()1969bl ). we derive the equa- 
tions of motion of the components of an eccentric mass- 
transferring binary with respect to a right handed inertial 
frame of reference OXYZ which has an arbitrary posi- 
tion and orientation in space (see Fig. [J) . We let Mi be 
the mass of star i at some time t at which mass is trans- 
ferred from the donor to the accretor, and Mi + SMi the 
mass of the same star at some time t + St, where St > 
is a small time interval. With these notations, SMi < 
corresponds to mass loss, and SMi > to mass accretion. 
We furthermore denote the point on the stellar surface 
at which mass is lost or accreted by Ai. For the donor 
star, Ai corresponds to the inner Lagrangian point Li, 
while for the accretor, Ai can be any point on the star's 
equator. For the remainder of the paper, we let « = 1 
correspond to the donor and i = 2 to the accretor. 

Because of the mass loss/gain, the center of mass of 
star i at time t + St is shifted from where it would have 
been had no mass transfer taken place. To describe this 
perturbation, we introduce an additional right-handed 
coordinate frame OiXiYiZi with a spatial velocity such 
that its origin follows the unperturbed orbit of star «, i.e., 
the origin of OiXiYiZi follows the path the center of mass 
of star i would have taken had no mass transfer occurred. 
Thus, at time t, the center of mass of star i lies at the 
origin of Oi, while at time t+ St it has a non-zero position 
vector with respect to Oi. We furthermore let the Zj-axis 
of the OiXiYiZi frame point in the direction of the or- 
bital angular momentum vector, and let the frame rotate 
synchronously with the unperturbed orbital angular ve- 
locity of the binary in the absence of mass transfer. The 
direction and orientation of the Xj-axes are then chosen 
such that at time t the A"i-axis points along the direction 
from the mass center of star i to the mass center of its 
companion. 

To describe the shift in the mass center of star i due to 
the mass loss/gain, we denote the position vector of Oi 
at times t and t + St with respect to the inertial frame 
by and R'i, respectively. The position vector of the 
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Fig. 1. — Schematic representation of the reference frames and position vectors adopted in the derivation of the equations of motion of 
the components of an eccentric mass-transferring binary. The Z- and Zi-axes of the OXYZ and OiXiYiZi frames are perpendicular to 
the plane of the page and are therefore not drawn. The geometry of the system at time t is shown in black, while the geometry at time t + St 
is shown in gray. The solid line connecting the origin of the 0\X\Y\Z\ frame at time t to the origin of the frame at time t + 5t represents 
the path the donor would have taken had no mass transfer occurred. The dashed line, on the other hand, represents the perturbed motion 
of the donor's mass center (small solid circles) due to the mass transfer. A similar perturbation is imparted to the motion of the accretor. 
For clarity, this perturbed motion and the O2X2Y2Z2 reference frame connected to the accretor (see text) are omitted from the figure. The 
dotted line illustrates a possible path of a mass element transferred from the donor to the accretor. The element leaves the donor at the 
inner Lagrangian point L\ (asterisk) and accretes onto the companion at the point A2 (open star). 



center of mass of star i at times t and t + 5t is then given 
by Ri and R[ + Jf/, where Jr/ is the position vector of 
the center of mass of star i at time t + 5t with respect 
to Oi. Moreover, we denote by rvi. and r^. the vectors 
from Oi to the position where the point Ai on the stellar 
surface would be at times t and t + 5t, respectively, had 
no mass been lost / accreted. The various position vectors 
at time t + 5t are related by 



(M, + 5 Ah) {r'^ +5rl]^ M, R- + SM, 



which, at the lowest order of approximation in SM, 
yields 

or, — ——— r 



(1) 
and 



5r' 



M,- 



As expected, the displacement of the center of mass of 
star i due to the mass loss/gain is directed along the line 
connecting the center of mass of the star and the mass 
ejection/accretion point. 



We furthermore denote with p^. the vector from the 
center of mass of star i at time t + 5t io the position 
where the point Ai would be at time t + 5t, had no mass 
been transferred between the binary components, and 
with 5p'j^, the perturbation of this vector caused by the 
mass transfer. It then follows that p'j^. ^ r^. — 5rl and 
thus, by definition, (5p^. = (5r/. At the lowest order of 
approximation in 5 Mi and <5r/, equation ([2]) therefore 
also yields 

, 5Mi ^, 
^PA.-^r^, (3) 

The definitions of and the relations between these various 
position vectors are illustrated schematically in Fig. [T] 

Next, we denote the absolute velocity of the center 
of mass of star i with respect to the inertial frame of 
reference at times t and t-\-5thy Vi and , respectively, 
and the absolute velocity of the ejected/accreted mass 

element by WsMi- The linear momentum Qi of star 1 at 
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time t is then given by 

Qi^MiVi, (4) 

and the total Hnear momentum Q'l of star 1 and the 
ejected mass element at time t + 6t hy 

Q[ = (Ml + SMi)v; - SMi WsM, ■ (5) 

Similarly, the total linear momentum Q2 of star 2 and 
the mass element to be accreted at time t is given by 

Q2= M2V2+SM^W5M^, (6) 

and the linear momentum Q'^ of star 2 at time t + 6t hy 

Q'2 = (M2 + 5M2)Vl (7) 

At time t + 5t the velocity of the center of mass of star 
i can be written as 



f^orb+'^fl^rb) Xfrt, (8) 



where Vk, is the absolute velocity of the origin of 



OiXiYiZi at time t + St, ^'^^^^ is the orbital angular ve- 
locity of the binary at time t + 6tm the absence of mass 
transfer, and Sfl'^^.^^ is the perturbation of the orbital an- 
gular velocity at time t + St due to the mass loss/gain of 
the binary components. 

In the limit of small 6t, taking the difference between 
the linear momenta and Qi, dividing the resulting 
equation by St, and noting that the absolute velocity 
of the origin of OiXiYiZi at time t is equal to Vi, yields 

M,^^F, + M,UsAu. (9) 
at 

Here, Fi = dQi j dt is the sum of all external forces acting 
on star i. Mi = dMi/dt is the mass loss/accretion rate 
of star i, and 

UsMi = WsAU - - dovh X rA, (10) 
is the relative velocity of the ejected/accreted mass el- 
ement with respect to the ejection/accretion point Ai. 
In the derivation of equation pO)) . we have made use of 
equations ^ and ^ and restricted ourselves to first- 
order terms in the small quantities SMi and Sfl'^^y^. 

The absolute acceleration of the center of mass of star 
i with respect to the inertial frame OXYZ is given by 

d^R 

= io, + 7rcl.i + 7cor,i, (H) 

where ~ dVOi I dt is the acceleration of the origin of 
OiXiY^Z, with respect to OXYZ, -frei,! = (Mj/M,)r^. 
is the relative acceleration of the center of mass of star i 
with respect to Oi, and 7cor,i = 2 {Mi/ Mi) (fiorb x rvij 
is the Coriolis acceleration of the center of mass of star i 
with respect to Oi^. The expressions for 7roi.i and 7cor,i 
follow from the observation that dp^./dt = (Mi /Mi) VAi, 
which one obtains by dividing equation ([3]) by St in the 
limiting case of small St. The equation of motion for the 
mass center of star i with respect to the inertial frame 
OXYZ then becomes 

d^R ^ • / ^ 

:Fi + Mi (UsM, + 2 (iorh X -TA, 



M,, 



dt^ 



(12) 



^ The centrifugal acceleration does not play a role since it is 
proportional to <5r.' which vanishes for small St. 



3.2. Relative Motion of the Binary Components 

We can now obtain the equation describing the rela- 
tive motion of the accretor (star 2) with respect to the 
donor (star 1) by taking the difference of the equations of 
motion of the stars with respect to the inertial frame of 
reference. For convenience, we first decompose the sum 
of the external forces acting on each star as 



^ GM^M2 Rr , - 

\r\ \Ri\ 



(13) 



where G is the Newtonian constant of gravitation, and 
fi the total gravitational force exerted on star i by the 
particles in the mass-transfer stream. It follows that 



d r 
di2 



G [Ml + M2) ^ /2 /i 
^3 '■+Af2 Ml 



^2 ,0 

■ T-J- I VSM2 + "orb X rA2 
M2 \ 

■ -JT V6M1 + "orb X VAi 
1V1\ V 

M2 ^ Ml ^ 



(14) 



M2 



Ml 



where r — R2 ^ Ri is the position vector of the accretor 
with respect to the donor, and vsMi = WsMi — ^ is 
the velocity of the ejected/accreted mass element with 
respect to the mass center of the mass losing/gaining 
star. 

Equation (jl4p can be written in the form of a perturbed 
two-body problem as 



G (Ml -I- M2) 



dt^ 



r + Sx + Ty + Wz, (15) 



where x is a unit vector in the direction of r, y is a 
unit vector in the orbital plane perpendicular to r in the 
direction of the orbital motion, and z is a unit vector 
perpendicular to the orbital plane parallel to and in the 
same direction as f^orb- The functions S, T, and W are 
found by taking the dot product of the perturbing force 
arising from the mass transfer between the binary com- 
ponents and the unit vector in the x, y, and z directions, 
respectively. These vector components are 



5 = 



f2,x fl.x 



Ml 

VSMi,x - 



M2 

Ml 
Ml 

h,y /: 



M2 
' M2 
M2. 
M2 
M 



VSAI2,. 



TA^ I cos ( 



If^orbllrAal Sini 

Ml , ^ 



Ml 



(16) 



^=M^-Mr + M2H-^ + "" 



rblr^al COS( 



^ (w5Afi,y + |^1orb||r^4l|) +^|rA2|sin0, (17) 

!l,z fl,z 



M2 



Ml ' 



(18) 



where (j) is the angle between x and the vector from the 
center of mass of the accretor to the mass accretion point 
A2, and the subscripts x, y, and z denote vector com- 
ponents in the x, y, and z directions, respectively. In 
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working out the vector products llorb x t^a^ , we assumed 
that A2 is located on the equator of the accreting star^. 
The terms contributing to the perturbed orbital motion 
can be categorized as follows: (i) term proportional to fi 
represent gravitational perturbation on the binary com- 
ponents caused by mass elements in the mass-transfer 
stream; (ii) terms proportional to Mi represent linear 
momentum exchange between the mass donor and accre- 
tor; and (iii) terms proportional to Mi represent shifts in 
the position of the mass centers of the mass donor and 
accretor due to the non-spherical symmetry of the mass 
loss or gain. In the limiting case where both stars are 
treated as point masses (Ir^iJ and Ir^j] — > 0), the 
only non-zero terms in the perturbed equations of motion 
are those due to gravitational perturbations of the mass 
transfer stream and the transport of linear momentum. 

3.3. Comparison with Previous Work 

The most recent study on the orbital evolution of ec- 
centric mass-transferring b i naries has been presented by 
iMatese fc Whitmird (iQSl flQSl hereafter MW83 and 
M W84, re s pectiv el y). These au thors extended the work 
of jHuand i|1956( l. iKruszewskil (jl964l ). and iPiotrowskil 
()1964[ ) by accounting for the effects of linear momen- 
tum transport between the binary components, as well 
as possible perturbations to the orbital motion caused 
by the mass transfer stream. However, they also derived 
the equations describing the motion of the binary com- 
ponents with respect to a frame of reference with origin 
at the mass center of the binary, which, for mass trans- 
ferring systems, is not an inertial frame of reference. The 
equations therefore do not account for the accelerations 
of the binary mass center caused by the mass transfer. 
Here, we demonstrate that if the procedure adopted by 
Matese & Whitmire is developed with respect to an in- 
ertial frame of reference that is not connected to the bi- 
nary, the resulting equations are in agreement with those 
derived in § 13.21 

The core of Matese & Whitmire's derivation is pre- 
sented in Section II of MW83. While the authors choose 
to adopt a reference frame with origin at the binary mass 
center early on in the investigation, the choice of the 
frame does not affect the derivation of the equations of 
motion up to and including their equation (24). In par- 
ticular, equation (13) in MW83, which, in our notation, 
reads 



= M,R, - M^fA, 



(19) 



is valid with respect to any inertial frame of reference 
with arbitrary position and orientation in space. The 
same applies to equation (1) in MW84: 



GM^M2 ,f' +h + ^^. (20) 



..|3 



In these equations, pi is the linear momentum of star i, 

and = Mi[Ri + vsai-) is the amount of linear mo- 
mentum transported by the transferred mass per unit 

^ For brevity, we refer to the point A2 as lying on the stellar 
surface. Though, in practice, it can lie at any point near the star 
where the transferred mass can be considered to be part of the 
accretor. For instance, if an accretion disk has formed around the 
accretor, it would be equally valid to write A2 as the point where 
the transferred mass impacts the outer edge of the accretion disk. 



time (see equation (3) of MW84). Substitution of equa- 
tion (I19p into equation ([20|) then yields 



M,R,,=-GM. 



Ri — R 



3-i 



3-i 



fi Mi / ■ \ M, ^ , , 



and thus 
lit? ' 



G(A/i+M2) A _ A 
1^1^ M2 Ml 

M2 1, \ Ml ( ^ 

+ j^^[v,M.^rA:,)- — [v,M.^rA, 



'2 

M2 ^ Ml ^ 



(22) 



where r — R2 — R\- Setting tvi . = iloib x , this equa- 
tion is in perfect agreement with equation (jl4p derived 
in §[321 

In MW83 and MW84, the authors incorrectly set 
^1 = -M2r/{Mi + M2) and R2 = Mir /{Mi + M2) 
in equation (j2ip . which is valid only when the origin of 
the frame of reference coincide with the mass center of 
the binary. For a mass-transferring binary, such a frame 
is, however, not an inertial frame and can therefore not 
be used for the derivation of the equations of motion of 
the binary components. Instead of equation (f^^ . Matese 
& Whitmire therefore find equations (7)-(8) in MW84, 
which lack the terms associated with the acceleration of 
the binary mass center due to the mass transfer. 

4. ORBITAL EVOLUTION EQUATIONS 

4.1. Secular Variation of the Orbital Elements 

In the classical framework of the theory of osculating 
elements, the equations governing the rate of change of 
the orbital semi-major axis a and eccentricity e due to 
mass transfer are obt ained from the perturbing functions 
S and T as (see, e.g.. 'Sternc"1960': 'B rouwer fc Clemencd 
ll96lHDanbv 1962; Fitzpatrick 1970) 



da 

Itt " n(l-e2)i/ 



-^[Sesinv + T{l + ecosv)], (23) 



dt na 
X I Ssinv + T 



2 cos I' + e (1 -I- cos^ v) 



1 



, (24) 



where n — 27r/Porb is the mean motion and v the true 
anomaly. These equations are independent of the per- 
turbing function W which solely appears in the equations 
governing the rates of change of the orbital inclination, 
the longitude of the ascending node, and the longitude 
of the periastron. 

After substitution of equations (fTB]) and (fTT]) for S and 
T into equations and (|24p . the equations governing 
the rates of change of the semi-major axis and eccentric- 
ity contain periodic as well as secular terms. Here we 
are mainly interested in the long-term secular evolution 
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of the orbit, and so we remove the periodic terms by 
averaging the equations over one orbital period: 



da 
'dt 

de\ 



1 



da ^ 
-^dt. 



-Poib ■/-P<„b/2 



— dt. 



Porb/2 



(25) 



(26) 



The integrals in these definitions are most conveniently 
computed in terms of the true anomaly, v. We therefore 
make a change of variables using 

dt = -4 — du. (27) 

n(l + ecosi^)^ 

For binaries with eccentric orbits, the resulting integrals 
can be calculated analytically only for very specific func- 
tional prescriptions of the mass-transfer rate Mi (e.g., 
when Ml is approximated by a Dirac delta function cen- 
tered on the periastron, see §[S|) . In general, the integrals 
must be computed numerically. 

4.2. Conservation of Orbital Angular Momentum 

Since the perturbing functions S and T depend on 
the properties of the mass transfer stream, calculation 
of the rates of secular change of the orbital semi-major 
axis and eccentricity, in principle, requires the calcula- 
tion of the trajectories of the particles in the stream 
(cf. iHadiidcmctriou 19693). long as no mass is lost 
from the system, such a calculation automatically incor- 
porates the conservation of total angular momentum in 
the system. Special cases of angular momentum conser- 
vation can, however, be used to bypass the calculation 
of detailed particle trajectories. Here, we adopt such a 
special case and assume that any orbital angular momen- 
tum carried by the particles in the mass-transfer stream 
is always immediately returned to the orbit, so that the 
orbital angular momentum of the binary is conserved. 

The orbital angular momentum of a binary with a 
semi-major axis a and eccentricity e is given by 

"Ga(l-e2)T^/' 



Jorb = M1M2 



Ml + M2 



so that 

<^orb 
«^orb 



Ml 
Ml 



M2 
M2 



IM1+M2 



1 a 
2a 



e e 



(28) 



(29) 



2 Ml -f- M2 2 a 1 - e2 

where a dot indicates the time derivative. 

In the case of eccentric orbits, substitution of equa- 
tions and into equation leads to 



^orb 
Jorh 



Ml 
M^ 



(1 



M2 1 Ml + M; 

M2 " 2M^^ 
- e2)i/2 



M2 



na{l + ecoav) 



■T. 



(30) 



As we shall see in the next section, by setting Mi + M2 = 

and < Jorb/>/orb >scc= and substituting equa- 
tion ([TT]) for T, equation ([50]) allows us to calculate the 
y-component of the final velocities of the accreting par- 
ticles as a function of their initial velocities without re- 
sorting to the computation of the ballistic trajectories of 
the mass transfer stream. 



In the limiting case of a circular orbit, equation (|29p 
is usually used to derive the rate of change of the semi- 
major axis of circular binaries under the assumption of 
conservation of both total mass (Mi = — M2) and orbital 
angular momentum (Jorb = 0): 



da Ml Ml 



(31) 



The assumption of orbital angular momentum conser- 
vation over secular timescales (< Jorb/ Jorb >scc= 0) 
is a standard assumption in nearly all investigations 
of conservative mass transfer in b inary systems (e.g., 
ISoberman et al.lfl997tlPribullalll998f l. which is valid over 
long timescales provided there is no significant storage of 
angular momentum in the spins of the components stars, 
the accretion flow, and/or the accretion disk. In future 
work, we will investigate the consequences of both mass 
and orbital angular momentum losses from the binary on 
the evolution of the orbital elements. 

5. ORBITAL EVOLUTION TIMESCALES 

In order to assess the timescales of orbital evolution 
due to mass transfer in eccentric binaries, we observe 
that, for eccentric binaries, mass transfer is expected to 
occur first at the periastron of the relative orbit, where 
the component stars are closest to each other. We there- 
fore explore the order of magnitude of the timescales as- 
suming a delta function mass transfer profile centered at 
the periastron of the binary orbit 



Ml = Mo 6 (ly) 



(32) 



where Mq < is the instantaneous mass transfer rate, 
and S{iy) is the Dirac delta function. 

We calculate the rates of secular change of the orbital 
semi-major axis and eccentricity from equations (j23p - 
((27|) . and neglect any gravitational attractions exerted 
by the particles in the mass-transfer stream on the com- 
ponent stars. Hence, we set 

/i,x = /2,. = 0, (33) 
/i,,-/2,,-0. (34) 

Substituting equations (HZ]), and ^ - ^ into 
equation ([30|) for < Jorb/ Jorb >soc= then yields a re- 
lationship between the initial and final y-component of 
the velocities of the transferred mass and the initial and 
final positions of the transferred mass given by 

qV5M2..y + '"5Muy = na(l - q) ( ) - \^orh,p\\fAi,P\ 



1 



-gl^^orb.pllrAsI C0S(/)p 1 



d<t) 

dv 



v=0 



(35) 



where the subscript P indicates quantities evaluated 
at the periastron of the binary orbit, q — M1/M2 is 
the binary mass ration, and we have used the relation 
dd{i>)/dv = —8(y)lv. Assuming the transferred mass 
elements are ejected by star 1 at the Li point with a ve- 
locity VfiMx equal to the star's rotational velocity at Li, 
we write 



•b,p||r4i,p|- 



(36) 
(37) 
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Moreover, under the assumption that each periastron 
passage of the binary components give rise to an ex- 
tremum of 4>{i'), the derivative dc/f/dvl^^Q is equal to zero 
in equation (|35p . so that 



VSAh 



■b.P 



(1 



q 



r^o cos 0p 



(38) 



For a binary with orbital period Porb — 1 day, eccen- 
tricity e = 0.2, and component masses A/i = 2 Mq and 
M2 — 1.44 M0, with |rAi.p| the distance from star 1 to 
Li (See Appendix A), and jrvijlcos^p « 2.9 x 10^ km 
(the circularization radius arou nd a compact object 
for th ese binary parameters; see iFrank. King, fc Raind 
(|2002l )). the accreting matter has a y- velocity component 
of the order of ~ —45 kms^^. 

After substitution of equation (|77|) and equations ([55)1 " 
([55]) , the integrals in equations ([^5]) and for the rates 
of secular change of the orbital semi-major axis and ec- 
centricity can be solved analytically to obtain 



da \ 



de\ 



a Mo 



1 



TT Ml (1 -e2)i/2 



qe- 



FA2 



(39) 



(1 



2tt 

\rA,,p\ 



Ml 



FAo 



+ 2(l-e)(<7-l) 



(40) 



We note that for a delta-function mass transfer rate given 
by equation ([32]) the x-component of the velocities vsMi 
and VSM2 does not enter into the derivation of theses 
equations due to the sin v term in equations (|23p and 
((24|) . Furthermore, in the limiting case of a circular orbit, 
equation ([39|) reduces to equation ([3T|) . provided that Afi 
in that equation is interpreted as the secular mean mass 
transfer rate (Mi)^^^ = Mo/(27r). In Appendix B, we 
present an alternative derivations to equations (|39|) and 
PO]) in the limiting case where the stars are treated as 
point masses. 

The rates of secular change of the semi-major axis and 
orbital eccentricity are thus linearly proportional to the 
magnitude of the mass transfer rate at periastron. Be- 
sides the obvious dependencies on a, e, q, and Mi, the 
rates also depend on the ratio of the donor's rotational 
angular velocity fli to the orbital angular velocity floih.p 
at periastron through the position of the Li point, f^^ . A 
fitting formula for the position of the Li point accurate to 
better than 4% over a wide range of g, e, and ili/iloib.p 
is given by equation (jA15|) in Appendix[Xl While the fit- 
ting formula can be used to obtain fully analytical rates 
of secular change of the semi-major axis and eccentric- 
ity, we here use the exact solutions for the position of Li 
obtained by numerically solving equation (|A13|) in Ap- 
pendix [Aj For a detailed discussion of the properties of 
the Li p oint in eccentric binari es, we refer the interested 
reader to ISepinskv et al.l (|2007D . 

To explore the effects of mass transfer on the orbital 
elements of eccentric binaries, we calculate the rates of 
secular change of the semi-major axis and eccentricity 
and determine the characteristic timescales = a/d and 



Te = e/e. While the actual timescales are given by the 
absolute values of and Tg, we here allow the timescales 
to be negative as well as positive in order to distinguish 
negative from positive rates of secular change of the or- 
bital elements. We also note that since |rvii,p| oc a (see 
Appendix A), the timescales do not explicitly depend on 
the orbital semi-major axis a except through the ratio 
Ir^^l/a of the radius of the accretor to the semi-major 
axis. For convenience, we therefore assume the accre- 
tor to be a compact object with radius jr^^l << a. The 
timescales are found to be insensitive to terms containing 
IrA^l/a in equation p9|) and (|40ll . Varying Ir^al from 
to O.Ola changes the timescales by less that 10%. In what 
follows, we therefore set Ir^al = 0. An implicit depen- 
dence on a may then still occur through the amplitude 
Mq of the mass transfer rate at periastron. Since incor- 
porating such a dependence in the analysis requires de- 
tailed modeling of the evolution of the donor star, which 
is beyond the scope of this investigation, we here restrict 
ourselves to exploring the timescales of orbital evolution 
for a constant Mq. The linear dependence of (d)^^^ and 

(^)scc on Mo in any case allows for any easy rescaling of 
our results to different mass transfer rates. 

In Fig. [21 we show the variations of Tq and Tg as 
functions of q for Mq — —10^^ Mqjt^^ and e = 
0.0, 0.1, . . . , 0.9. In all cases, the donor is assumed to 
rotate synchronously with the orbital angular velocity at 
the periastron, and the accretor is assumed to be a neu- 
tron star of mass M2 — 1.44 Mq. The timescales of the 
secular evolution of the semi-major axis show a strong 
dependence on q, and a milder dependence on e, unless 
e > 0.7. The timescales for the secular evolution of the 
orbital eccentricity always depend strongly on both q and 
e. These timescales can furthermore be positive as well 
as negative, so that the semi-major axis and eccentricity 
can increase as well as decrease under the influence of 
mass transfer at the periastron of the binary orbit. 

From Fig. [21 as well as equations (|39| and (j40|) . it can 
be seen that, for a given ratio of the donor's rotational 
angular velocity to the orbital angular velocity at peri- 
astron, the line dividing positive from negative rates of 
secular change of the orbital elements is a function of 
q and e. This is illustrated further in Fig. [3] where the 
timescales of orbital evolution are displayed as contour 
plots in the (g, e)-plane. The thick black line near the 
center of the plots marks the transition values of q and 
e where the rates of secular change of a and e transition 
from being positive (to the left of the thick black line) to 
negative (to the right of the thick black line). Varying 
^^i/^^oib,p between 0.5 and 1.5 changes the position of 
the transition line by less than 10% in comparison to the 
f^i/f^orb.p = 1 case displayed in Figures [21 and [3l 

In the limiting case of a circular orbit, the orbit ex- 
pands when g < 1 and shrinks when g > 1, in agreement 
with the classical result obtained from equation (|3T|) . For 
non-zero eccentricities, the critical mass ratio separat- 
ing positive from negative values of (d)^^^ decreases with 
increasing orbital eccentricities. This behavior can be 
understood by substituting the fitting formula for the 
position of the Li point given by equation (jAlsp in Ap- 
pendixIXIinto equation (|59|) and setting (d)^^,^ = 0. How- 
ever, we can fit the critical mass ratio separating expand- 
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Fig. 2. — Orbital evolution timescales for a delta function mass transfer profile centered at the periastron of the binary orbit with an 
instantaneous mass transfer rate of Mq = — lO~^M0yr~^. The timescales are calculated under the assumption that the donor rotates 
synchronously with the orbital angular velocity at periastron, and that the accretor is a 1.44 Mq neutron star. Shown at left (right) are the 
timescales for the evolution of the semi-major axis (eccentricity) as a function of the mass ratio, q, for a range of eccentricities, e. Regimes 
where the timescale is negative correspond to a decrease of the semi- major axis (eccentricity), while regimes where the timoscale is positive 
correspond to an increase of the semi-major axis (eccentricity). 
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Fig. 3. — Contour plot of the orbital evolution timescales for the semi-major axis and eccentricity in the (g, e)-plane for the same set of 
assumptions as adopted in figure[2] Timescales for the evolution of the semi-major axis are shown on the left; timescales for the evolution 
of the orbital eccentricity on the right. The diff'erent shades of gray designate regions of the (g, e) parameter space with timescales (in Gyr) 
in the ranges labeled in the plots. The thick black line near the center of each plot designates the transition point where the rate of change 
of the semi-major axis or orbital eccentricity changes from positive (to the loft of the thick black line) to negative (to the right of the thick 
black line). 



ing from shrinking orbits with a simpler formula given by 



9crit ^ 1 - 0.4e + 0.18e^ 



(41) 



The critical mass ratio separating positive from neg- 
ative values of (e)^^,^. is largely independent of e. Pro- 
ceeding in a similar fashion as for the derivation of equa- 
tion (j4ip . we derive the critical mass ratio separating 
increasing from decreasing eccentricities to be approxi- 
mately given by 



Last, we note that a more quantitative numerical com- 
parison between the above approximation formulae for 
(/crit and the exact numerical solutions shows that equa- 
tions and ([1^ are accurate to better than 1%. 



gcrit ^ 0.76 + 0.012e. 



(42) 



6. TIDAL EVOLUTION TIMESCALES 

A crucial question for assessing the relevance of the 
work presented here is how the derived orbital evolution 
timescales compare to the corresponding timescales as- 
sociated with other orbital evolution mechanisms such as 
tides. In Fig.[4l we show the secular evolution timescales 
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of the semi-major axis and orbital eccentricity of a mass- 
transferring binary due to tidal dissipation in the donor 
star as a function of q, for different values of the eccentric- 
ity, e. The timescale s are strong function s of jr^J/aand 
are determined as in iHurlev et al l (l20nl )3 ( 

see also lZahnI 

[T977l[l97a[H^[l98lh . The radius | r^i I is determined by 
assuming the donor is on the zero-age main sequence and 
that the orbital separation is th en obtained by eq uating 
the radius of the donor (given bv lTout et al.|[l996l ) to the 
volume-equivalent radius of its Roche lobe at the peri- 
astron of the binary orbit (see iSepinsky et al . 2007j) . As 
before, we assume the donor rotates synchronously with 
the orbital motion at periastron and that the accretor is 
a 1.44 M© neutron star. 

The timescales of orbital evolution due to tides range 
from a few Myr to more than a Hubble time, depending 
on the binary mass ratio and the orbital eccentricity. The 
discontinuity in the timescales at g ~ 0.87 corresponds 
to the transition from donor stars with convective en- 
velopes (Ml < 1.25 Mq) to donor stars with radiative 
envelopes (Mi > 1.25 M©) which are subject to different 
tidal dissipation mechanisms. It follows that tides do not 
necessarily lead to rapid circularization during the early 
stages of mass transfer, especially for orbital eccentrici- 
ties e > 0.3. Furthermore, for the adopted system pa- 
rameters, the orbital eccentricity always decreases, while 
the orbital semi-major axis can either increase or de- 
crease. Hence, in some regions of the parameter space, 
the effects of tides and mass transfer are additive, while 
in other regions they are competitive. This is illustrated 
in more detail in figure [5] where we show the orbital evo- 
lution timescales due to the combined effect of tides and 
mass transfer. In the calculations of the timescales, we 
have assumed that, at the lowest order of approximation, 
the effects of tides and mass transfer are decoupled. The 
total rate of change of the orbital elements is then given 
by the sum of the rate of change of the orbital elements 
due to tides and mass transfer. 

When q < 0.87 and e > 0.4, the effects of tides and 
mass transfer on the orbital semi-major axis are always 
opposed, with the orbital expansion due to mass trans- 
fer dominating the orbital shrinkage due to tides. In the 
case of the orbital eccentricity, the increase of the ec- 
centricity due to mass transfer dominates the decrease 
due to tides for mass ratios smaller than some critical 
mass ratio which depends strongly on the orbital eccen- 
tricity. Since the timescales of orbital evolution due to 
mass transfer are inversely proportional to the magnitude 
A/q of the mass-transfer rate at periastron, the param- 
eter space where and the extent to which mass transfer 
dominates increases with the rate of mass transfer at pe- 
riastron. 

In figure [H we show the total orbital evolution 
timescale due to the sum of tidal and mass transfer ef- 
fects as a contour plot in the (9, e)-plane. The thick 
black lines indicate the transitions from positive (left 
of the thick black line) to negative (right of the thick 

^ N ote that there is a typo in equation (42) of IHurlev et al.l 
l|2002l ) . The correct equation for k/T for stars with radiative en- 
velopes is (J. Hurley, Private Communication) 

(fc/T)^ = 1.9782 X 10^ {MR''/a''f'^ (1 + q-^f^ E2 yr"!. 



black line) rates of change of the semi-major axis and 
eccentricity. The white dividing line near q « 0.87 corre- 
sponds to the transition between tidal dissipation mecha- 
nisms in stars with convective envelopes (Mi < 1.25 M©) 
and stars with radiative envelopes {Mi > 1.25 M©). It 
follows that there are large regions of parameter space 
where the combined effects of mass transfer and tidal 
evolution do not rapidly circularize the orbit. In partic- 
ular, for q > 0.87 and e > 0.75 orbital circularization 
always takes longer than 10 Gyr, while for q to the left of 
the thick black line the orbital eccentricity grows rather 
than shrinks. For a given q left of the thick black line, 
the timescales for eccentricity growth increase with in- 
creasing e though, so that there is no runaway eccen- 
tricity growth. Hence, for small q, mass transfer at the 
periastron of eccentric orbits may provide a means for 
inducing non-negligible eccentricities in low-mass binary 
or planetary systems. The orbital semi-major axis, on 
the other hand, always increases when e < 0.55, but can 
increase as well as decrease when e > 0.55, depending on 
the binary mass ratio q. We recall that both the tidal 
and mass transfer orbital evolution time scales depend 
on the ratio of the donor's rotational angular velocity fii 
to the orbital angular velocity f^orb.p and that we have 
set r2i/r2orb,p = 1 in all figures shown. 

7. CONCLUDING REMARKS 

We developed a formalism to calculate the evolution 
of the semi-major axis and orbital eccentricity due to 
mass transfer in eccentric binaries, assuming conserva- 
tion of total system mass and orbital angular momentum. 
Adopting a delta-function mass-transfer profile centered 
at the periastron of the binary orbit yields rates of secu- 
lar change of the orbital elements that are linearly pro- 
portional to the magnitude Mq of the mass-transfer rate 
at the periastron. For Mq = 10~^ Mq yr"-'^, this yields 
timescales of orbital evolution ranging from a few Myr 
to a Hubble time or longer. Depending on the initial 
binary mass ratio and orbital eccentricity, the rates of 
secular change of the orbital semi-major axis and eccen- 
tricity can be positive as well as negative, so these orbital 
elements can increase as well as decrease with time. 

Comparison of the timescales of orbital evolution due 
to mass transfer with the timescales of orbital evolution 
due to tidal dissipation shows that the effects can either 
be additive or competitive, depending on the binary mass 
ratio, the orbital eccentricity, and the magnitude of the 
mass-transfer rate at the periastron. Contrary to what 
is often assumed in even the most state-of-the-art binary 
evolution and population synthesis codes, tides do not al- 
ways lead to rapid circularization during the early stages 
of mass transfer. Thus, phases of episodic mass transfer 
may occur at successive periastron passages and may per- 
sist for long periods of time. As a first approximation, the 
evolution of the orbital semi-major axis and eccentricity 
due to mass transfer in eccentric binaries can be incor- 
porated into binary evolution and population synthesis 
codes by means of equations ([39|) and (|40|) in which the 
mass-transfer rate is approximated by a delta-function 
of amplitude Mq centered at the periastron of the binary 
orbit. 

In future papers, we will relax the assumption of con- 
servation of total system mass and orbital angular mo- 
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Fig. 4. — Timescales of orbital evolution due to tidal dissipation in a Roche-lobe filling component of a close binary under the assumption 
that the donor is a zero-age main-sequence star rotating synchronously with the orbital angular velocity at the periastron, and the accretor 
is a 1.44 Mq neutron star. Shown at left (right) are the timescales for the evolution of the semi-major axis (eccentricity) as a function of 
the mass ratio, q, for different orbital eccentricities e. Regimes where the timescales are negative correspond to a decrease of the semi-major 
axis (eccentricity), while regimes where the timescales are positive correspond to an increase of the semi-major axis (eccentricity). The 
discontinuity at q ~ 0.87 corresponds to the transition from donor stars with convective envelopes (Mi < 1.25 M©) to donor stars with 
radiative envelopes (Mi > 1.25 M©). 
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Fig. 5. — Orbital evolution timescales due to the combined effects of tidal dissipation in a Roche Lobe filling component of a close 
binary system and a delta function mass transfer with an amplitude Mq = —10^'* Af© yr^"^ centered at the periastron of the orbit under 
the assumption that the donor is rotating synchronously with the orbital angular velocity at periastron, and the accretor is a 1.44 M© 
neutron star. The contribution to the orbital evolution timescales due to tides is determined under the assumption that the donor is 
zero-age main-sequence star. Shown at left (right) are the timescales for the evolution of the semi-major axis (eccentricity) as a function 
of the mass ratio, q, for a range of eccentricities, e. Regimes where the timescale is negative correspond to a decrease in the semi-major 
axis (eccentricity), while regimes where the timescales is positive correspond to an increase in the semi-major axis (eccentricity). The 
discontinuity at y ~ 0.87 corresponds to a transition of the dominant tidal dissipation mechanism which is different for donor stars with 
convective envelopes (Mi < 1.25 Mq) than for donor stars with radiative envelopes (Mi > 1.25 M©). 
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Fig. 6. — Contour plots of the total orbital evolution timescales in the (g, e)-plane due to the combined effects of tidal dissipation and 
mass transfer for the same set of assumptions adopted in figure [5] Timescales for the evolution of the semi-major axis are shown at 
left; timescales for the evolution of the orbital eccentricity are shown at right. The different shades of gray designate regions of the (g, e) 
parameter space with timescales (in Gyr) in the ranges labeled in the plots. The thick black line designates the transition point where the 
rate of change of the semi-major axis or orbital eccentricity changes from positive (to the left of the thick black line) to negative (to the 
right of the thick black line). The vertical white line at g ~ 0.87 corresponds to a transition of the dominant tidal dissipation mechanism 
which is different for donor stars with convective envelopes (A/i < 1.25 Mq) than for donor stars with radiative envelopes (A/i > 1.25 AIq). 
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Fig. 7. — Schematic diagram showing the vectors pertinent to the derivation of the equipotential surfaces for a non-synchronous, eccentric 
binary. The center of mass of the system is shown as a cross near the center of the diagram and the center of mass of each star is shown 
as a filled circle at the star's center. The vectors _Ri and R2 begin at the center of mass of the system. 



mentum, and examine the effects of non-conservative 
mass transfer on the orbital elements of eccentric bina- 
ries. We also intend to study the onset of mass transfer in 
eccentric binaries in more detail, adopting realistic mass- 
transfer rates appropriate for atmospheric Roche-lobe 



overflow in interacting binaries as discussed by iRitteil 
(Il988). We will consider individual binary systems that 
are known to be eccentric and transferring mass during 
periastron passage, as well as populations of eccentric 
mass-transferring binaries and their descendants. 



APPENDIX 

EQUIPOTENTIAL SURFACES AND THE INNER LAGRANGIAN POINT IN ECCENTRIC BINARIES 

A crucial element in the description of mass transfer in any binary system is the location of the inner Lagrangian 
point (Li) through which matter flows from the donor to the accretor. While a solution for the location of Li is not 
analytic, the case of circular orbits with synchronized components can be approximated by the formula 

= 0.5 + 0.22 log g, (Al) 

where X^-^^ is the distance of the Li point from the mass center of the donor star in units of the distance between 
the stars, and q — M1/M2 is the mass ratio of the binary defined as the ratio of the donor mass Mi to the accretor 
mass M2. For mass ratios in t he range 0.1 < g < 15, this formula is valid to an accuracy of bet t er tha n 2% (e.g. 
iDrobvshevski fc Reznikovlll974f ). The above formula has been generalized bv lPratt fc Strittmatterl (|1976( ) to include 
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the effect of a non-synchronously rotating donor star: 



X 



Li 



0.53-0.03 



orb 



l + ^logg). (A2) 



Here, fli is the rotational angular velocity of the donor, and fiorb the orbital angular velocity of the binary. For 
< r^i/Oprb < 1 and 0.2 < q < 10, the formula is vahd to an accuracy of better than 3%. It is to be noted though that 
iPratt fc Strittmatteil (|1976f ) derived the position of the Li point under the assumption that the orbital and rotational 
periods of the binary and its component stars are much longer than the dynamical timescale of the donor. Despite 
the better than 3% accuracy of the fit given by equation (jA2[) . the formula may therefore still break down because the 
underlying formalism is no longer valid. 

In this paper, it is necessary to generalize the formula for the position of the Li point even further to account for a 
non-zero eccentricity as well as non-synchronous rotation. For this purpose, we determine the equipotential surfaces 
describing the shape of the c o mpon ents of a non-synchronous, eccentric binary. Our procedure is a generalization of 
the steps outlined by lLimb"erl ()1963[ ) for the derivation of the equipotentials of a non-synchronous, circular binary. 

Here, we consider an eccentric binary system where the stars are considered to be centrally condensed and spherically 
symmetric, and thus can be well described by Roche models of masses Mi and M2. Their orbit is assumed to be 
Keplerian with semi-major axis a and eccentricity e. Star 1 is furthermore assumed to rotate uniformly with a 
constant angular velocity ili parallel to the orbital angular velocity floib- We note that for an unperturbed Keplerian 
orbit, the magnitude of fiorb is a function of time, while the direction and orientation remain fixed in space. 

As is customary, we determine the equipotential surfaces with respect to a Cartesian coordinate frame OXYZ with 
origin O at the mass center of star 1, and with the Z-axis pointing along The X- and F-axes are co-rotating 
with the star at angular velocity Qi. Since non-synchronously rotating binary components are inevitably subjected 
to time-dependent tides invoked by their companion, the mass elements will oscillate with frequencies determined by 
the difference between the rotational and orbital frequencies. Here, we neglect these tidally induced oscillations as 
well as any other type of bulk motion of matter due to, e.g., convection. This approximation is valid as long as the 
characteristic timescale associated with the motion of the mass elements is sufficiently long compar ed to the star 's 
dynamic al timescale (more extended discussions on the validity of the approximation can be found in iLimbeil (|1963[ ). 
ISavonii j ((19781 . and lSepinskv et ahl ([2007)). Star 1 is then completely stationary with respect to the co-rotating frame 
of reference. 

With respect to the OXYZ frame, the equation of motion of a mass element at the surface of star 1 is given by 



ri = fo - i?i - Oi X (^Oi X fi j - 2 fii X ri, (A3) 

where tq and fi are the position vectors of the mass element with respect to the mass center of the binary and star 
1, respectively, and Ri is the position vector of the center of mass of star 1 with respect to the center of mass of the 
binary, and — ili x (f2i x fi) and — 2ili x fi are the centrifugal and Coriolis acceleration with respect to the rotating 
coordinate frame. 

Under the assumption that, in an inertial frame of reference, the only forces acting on the considered mass element 
are those resulting from the pressure gradients in star 1 and the gravitational attractions of star 1 and its companion, 
the acceleration of the mass element with respect to the binary's center of mass is given by 

^o = --VP-vf-G^-G^), (A4) 
P \ \r2\J 

where G is the Newtonian constant of gravitation, p the mass density, P the pressure, and r2 the position vector of 
the mass element with respect to the mass center of star 2. The gradient in equation (jA4p and subsequent equations 
in this Appendix are taken with respect to X , Y, and Z. 

The acceleration of the center of mass of star 1 with respect to the binary mass center, on the other hand, is given 
by 

A GM2 Ri , . , 

where D{t) (simplified to D in what follows) is the time-dependent distance between star 1 and star 2. 
Substituting equations (|A4[) and (jASp into equation (jA3D . we obtain 

h = --vp-v(-G^-G^) + ^^-nix{nixn)^2nixh. (A6) 

When the X-axis coincides with the line connecting the mass centers of the two stars, the third and fourth term in 
the right-hand member of equation (jA6p can be written as the gradient of a potential function as 

GM,Ri_ _ . (GM2 \ 
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(A8) 



where X and y are the Cartesian coordinates of the mass element under consideration. With these transformations, 
equation (jA6p becomes 

?i = --VP -VVi -2ni xr'i, (A9) 
P 



V, = -G^, - Gp-^ \m\X^ + Y') + ^X. (AlO) 



where 



Since star 1 is assumed to be static in the rotating frame, fi^fi — 0, so that equation (|A9p reduces to 

VP = -pVVi. (All) 

This equation governs the instantaneous shape of the equipotential surfaces at the instant when the X-axis coincides 
with the hne connecting the mass centers of the two stars. However, since this instant is entirely arbitrary, the equation 
is generally valid and can be used to determine the shape of the star at an y phase of th e or bit, with appropriate re- 
definition of the X-axis at every instant. A similar equation was derived bv lAvnil (|1976l ) and lWilsonI (jl979). 

Next, we look for the critical equipotential surface for which matter starts to flow from star 1 to star 2 through 
the inner Lagrangian point Li. For circular binaries with synchronous ly rotating component star s, the Li point is 
always located on the line connecting the mass centers of the two stars. iMatese fc Whitmird ()1983l ) have shown that 
non-synchronous rotation combined with a misalignment between the spin axis of star 1 and the orbital plane may 
cause the Li point to oscillate in the Z-direction. As a first approximation, we here neglect these oscillations and 
assume the Li point to be located on the line connecting the mass centers of the two stars. The position of the Li 
point is then obtained by setting dVi/dX = 0. The resulting equation for the position of Li is 

which can be rewritten in dimensionless form as 

^ -/^XL,(l + g)(l + e) + l = 0. (A13) 

Here, X^-^ — X/D is the position of the Li point on the X-axis in units of the distance between the two stars, 
q = M1/M2 is the binary mass ratio, / = Sli/ilorb.p is the ratio of the star I's rotational angular velocity to the 
orbital angular velocity at periastron, and ^^^i, p = G{Mi + M2)(l -I- e)/Dp. In the latter expression. Dp = a{l — e) 

is the distance between the stars at periastron. Note that in th e particular case o f a c ircular orbit, equation (IA13|1 
reduces to the equation for the position of the Li point derived bv lKruszewskil ()1963f ) and lPratt fc Strittmatterl (|l976). 

Equation (jAlSp must be solved numerically for the location of Li. Since we have made no explicit assumptions 
about the relative location of the two star in the orbit, and we have let the distance between the stars D ~ a{l ~ 
e^)/(l + ecosv), with v the true anomaly, be an explicit function of time, this equation characterizes the position of 
the Li point at any position in the orbit. Figure [8] shows the relative change 

XlA'^^Q) ^ ' 

in the position of the Li point as a function of the true anomaly for g = 1, / = 1, and a variety of eccentricities. 
The position of the Li point over the course of a single orbit varies by less than 15% when e < 0.1, and by 30-40% 
when e > 0.2. At apastron, Axj^i increases with increasing eccentricity when e < 0.5, and decreases with increasing 
eccentricity when e > 0.5. This turn-around at e ~ 0.5 is caused by the increasing deviations from synchronous 
rotation at apastron with increasing orbital eccentricity. 

We have also derived a fitting formula for the position of the Li point at the periastron of the binary orbit as a 
function of q, e, and /: 

Xl, = 0.529 + 0.231 logg-/2(0.03H-0.025e)(H-0.41og(7). (A15) 

This fitting formula is accurate to better than 4% for mass ratios 0.1 < <? < 10, eccentricities < e < 0.9, and ratios 
of stellar rotational angular velocities to orbital angular velocities at p eriastron < / < 1.5. In the limiting case of a 
circular, synchronized binary, equati on ()A15|) reduces to equation (jAip . while for a circular, non-synchronized binary, 
equation (|A15p reduces to equation (|A2[1 . A detailed discussion of the properties of aU five L agrangian points and the 
volume of the critical Roche Lobe in eccentric binaries is presented in ISepinskv et al.l (|2007[ ) 
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Fig. 8. — The relative change in the position of the Li point as a function of true anomaly u, for g = 1, / = 1, and a range of eccentricities. 



ORBITAL EVOLUTION DUE TO INSTANTANEOUS MASS TRANSFER BETWEEN TWO POINT MASSES 

When the components of a mass-transferring binary are treated as point masses and mass transfer is assumed to 
be instantaneous, an alternative derivation of the equations governing the secular change of the orbital semi-major 
axis and eccentricity may be obtained from the equations for the specific orbital energy and angular momentum. This 
method is commonly used to relate the post- to pre-supernova orbital parameters of binaries in which one of the 
comp onents undergoes an instantaneous supernova explosion (e.g. lHillslll983t iBrandt fc Podsiadlowski|[l995l : iKalogeral 
Il996l ). As in Section [5l we assume conservation of total system mass and describe the mass transfer by a delta function 
mass-transfer rate centered at the periastron of the binary orbit [see equation (j32|) ]. 

The relative orbital velocity, Koi, of two stars with masses Mi (the donor star) and M2 (the accreting star) in a 
binary with orbital semi-major axis a and eccentricity e is given by 

V^,, = G(Mi + M2) (|j - ^) ' (Bl) 

where G is the Newtonian constant of gravitation, and f is the relative position vector of the accretor with respect 
to the donor. Under the assumption that mass elements are transferred instantaneously between the two stars, the 
distance |r| is the same right before and right after mass transfer, so that the rate of change of the orbital semi-major 
axis due to instantaneous mass transfer is given by 

da _ 2a^l/rei dV^d . 
dt ~ G{Mi +M2) dt ' ^ ' 
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Next, the specific orbital angular momentum of the binary written in terms of the binary component masses and 
orbital elements is given by 

|f X K-eip = G{Mi + A/2)a(l - e^). (B3) 



Noting that |r x V,ci| = |r||Vrci| at periastron, and proceeding in a similar way as for the derivation of equation (|B2p . 
we derive the rate of change of the orbital eccentricity due to instantaneous mass transfer to be 

de ^ 2a{l - e)\Vrci\ d\Vr,i\ 

dt G(Mi + Ah) dt ' ^ ' 



Finally, we eliminate Vtci from equations (jB2p and (jB4p to obtain equations for the rates of changes of the orbital 
semi-major axis and eccentricity in terms of the binary orbital elements and component masses. For this purpose, we 
use the conservation of linear momentum to write the rate of change of the relative orbital velocity due to mass-transfer 
as 

dtcx A/2 . Ml ^ 

^^M^'"'^"^"'''^' ^^^^ 
where v^Mi and wsa/s are the relative velocities of the transferred mass elements with respect to the mass center of the 
donor and accretor, respectively. To obtain an equation for the change in the magnitude of Vrei, we write the vector 
components of Vrd, vsm^, and vsm2 with respect to the x, y, and z unit vectors introduced in Section [321 and take 
the dot product of equation (|B5P with V^ei- Since, l^ei = V^ci y at the periastron of the binary orbit, it follows that 

~dr " "Ml '^"'5^^=''' + '"^'^-'^■y) ■ 
Assuming conservation of orbital angular momentum it then follows from equations (j29|) and (jB2p - (|B6P that 

1/2 



where we have used Kepler's 3rd law and that |K-ei| = (G(Mi + M2)/a)^^^ ((1 -I- e)/(l — e))^^^ at the periastron of 
the binary orbit. This equation is identical to equation ([35]) in the limiting case where |r/ij — *■ and Ir^ial ^ 0. 
Substituting equations (eq-C-vdm2) and (eq-C-vrel) into equations (|B2p and (|B4p and averaging over the orbital 
period, we find 

/ sec ^ 
/ sec ^ 

These equations are identical to equations ([25]) and ([26)1 in the limiting case where |rvij and |rA2l ~^ 0- 
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